FAR ULTRAVIOLET IMAGERY OF THE EDGE-ON SPIRAL 

GALAXY NGC 4631 
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ABSTRACT 

Far ultraviolet (FUV) imagery of the edge-on, Sc/SBd galaxy, NGC 4631 reveals very strong 
FUV emission, resulting from active star formation, uniformly distributed along the galactic mid- 
plane. Multi-band imagery, H I and H II position-velocity curves and extinction considerations 
all imply that the emission is from the outer edges of the visible galaxy. The overall FUV 
morphology of this edge-on disk system is remarkably similar to those of the so-called "chain 
galaxies" evident at high redshift, thus suggesting a similar interpretation for at least some of 
those distant objects. FUV, U, B and V magnitudes, measured for 48 star forming regions, 
along with corresponding Ha and H(3 measurements are used to construct diagnostic color-color 
diagrams. Although there are significant exceptions, most of the star forming regions are less 
massive and older than 30 Doradus. Comparison with the expectations from two star formation 
models yields ages of 2.7 to 10 Myr for the instantaneous burst (IB) model and star formation 
cut-off ages of to 9 Myr for the continuous star formation (CSF) model. Interpreted in terms 
of the IB model the photometry implies a total created mass in the 48 star forming regions of 
2.5x10^ M0. When viewed as resulting from constant star formation the photometry implies 
a star formation rate of 0.33 Mq yr~^. These results are compared to those derived from FIR 
and radio observations. Corrections for FUV emission reprocessed by interstellar grains are 
estimated. 

A large ring, ~3 kpc in diameter, of 14 star forming regions is concentrically located with 
an expanding H I shell (Rand & van der Hulst 1993) toward the eastern end of the galaxy. Our 
observations imply that the shell may have been generated primarily by supernovae arising from 
5.3x10* OB stars in a massive star forming region beginning about 20 Myr ago, and that the 
presently observed FUV bright emission is due to second generation stars. 



Subject headings: galaxies: individual (NGC 4631)- 
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1. Introduction 

The group of galaxies including NGC 4631 pro- 
vides an outstanding example of a galaxy interac- 
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tion accompanied by intensive star formation. Far 
Ultraviolet (FUV) imagery of the Sc/SBd galaxy 
NGC 4631, recorded by the Ultraviolet Imaging 
Telescope (UIT), exhibits very bright FUV emis- 
sion from sources corresponding to the H II re- 
gions cataloged by Crillon & Monnet (1969) along 
with substantial diffuse FUV light. The observed 
emission is, in our experience, extraordinary in 
that NGC 4631 is observed nearly edge-on, and 
strong attenuation of FUV light by dust in the 
galactic disk would be anticipated. Such is the 
case for NGC 891, another highly inclined galaxy 
and M82, an edge-on starburst galaxy - both 
of which show very little detectable FUV emis- 
sion. (Images of NGC 891 and M82 are accessible 
at http://archive.stsci.edu/astro/uit/index.html). 
Together with the Ha imagery obtained by Rand 
et al. (1992) the FUV results signify widespread 
intensive star formation along the observable disk 
of NGC 4631. The overall FUV morphology is re- 
markably similar to those of the so-called "chain 
galaxies" that have been observed at high redshift 
(Cowie et al. 1995). It may be that a significant 
fraction of these distant FUV sources is composed 
of highly inclined disk galaxies rather than dynam- 
ically unstable chains of sub-galactic clumps. See 
Smith et al. (1997) for further details. 

Observations of the radio halo imply that the 
galactic-wide star formation in NGC 4631 is suf- 
ficiently intense as to supply cosmic rays to the 
halo through bursting supershells formed by su- 
pernovae and stellar winds in the disk (Dahlcm et 
al. 1995; GoUa & Hummel 1994). Optical images 
reveal a patchy distribution of dust rather than a 
well defined dust lane. Thus, the extensive FUV 
emission may result from the combined dispersal 
of dust by these bursting supershells and by forces 
generated through tidal interactions NGC 4631 
and NGC 4656 (Combes 1978). 

There is also evidence for an inner disk or 
torus (r^50") revealed most clearly in CO(l-O) 
and C0(2-l) line emissions (GoUa & Wielebinski 
1994; Sofue et al. 1989, 1990). The radio con- 
tinuum emission coincides spatially with the CO 
emission and has a substantial non-thermal com- 
ponent (Duric et al. 1982). There also appears 
to be a region of enhanced star formation near 
the center (R<360 pc) of NGC 4631 (Duric et al. 
1982; Sofue et al. 1989) coinciding closely with the 
center of the galaxy as measured in the infrared 



by Aaronson (1978). More recently, Wang et al. 
(1995) reported the ROSAT detection of soft X- 
rays extending out to 8 kpc above the galactic mid- 
plane. When interpreted in terms of the burst su- 
pershell model, this observation implies substan- 
tial star formation near the center of NGC 4631. 
Indeed, considering mostly the radio data, Golla 
& Wielebinski (1994) view NGC 4631 as a galaxy 
with a mild central starburst. 

Of interest is a supershell with a diameter of 
3 kpc detected at 21 cm by Rand & van dcr Hulst 
(1993). The shell corresponds to a distortion in 
the disk near the eastern end of the galaxy appear- 
ing in all optical images and in a particularly obvi- 
ous way in both the FUV and Ha images. A pre- 
liminary explanation for the shell is that it is pro- 
duced by the combined effects of stellar winds and 
supernovae from OB associations containing (1- 
3.5) X 10'* OB stars (Rand & van der Hulst 1993). 
Alternatively, these authors suggest that the east- 
ern cavity may have been punched into the galactic 
disk by a high velocity cloud (HVC) which itself 
originated from the tidal stripping of the galactic 
disk during a prior interaction with NGC 4656. A 
collision with a dwarf galaxy could have also pro- 
duced the same result. 

In this paper we present the FUV imagery 
of NGC 4631. Using these data together with 
recently obtained optical imagery, we infer the 
star formation history of the bright FUV emit- 
ting associations including those circumscribed 
by the Rand and van der Hulst super shell. 
The implications of these results on the ultra- 
violet emission morphology, local and global 
star formation, and the viability of the star- 
burst shell and HVC impact hypothesis are dis- 
cussed. In all subsequent discussions, we place 
the galactic center at RA(2000) = 12'*39" 39?8, 
DEC(2000) = -f32°48'48", adopt a galactic dis- 
tance of 7.5 Mpc (Golla & Hummel 1994) and 
abbreviate the H II regions observed by Crillon & 
Monnet (1969) as CM(#). 

2. Observations 

NGC 4631 was observed by the UIT through 
the Bl filter (Ae// = 152.0 nm; AA = 35.4 nm) 
for 1140 seconds on March 10, 1995 during the 
Spacelab ASTRO-2 mission. The instrument is 
completely described in Stecher et al. (1992, 1997). 



FUV images, intensified and down-converted to 
visible wavelengths by a two stage magnetically 
focussed intensifier, were recorded on IIa-0 film, 
digitized with a PDS microdensitometer, and re- 
duced as described in Stecher et al. (1992, 1997). 
The plate scale of 1' mm~^ together with a micro- 
densitometer aperture of 0.021 mm resulted in a 
pixel size of 1.26". The FWHM of the point spread 
function was 3.4" over a circular 40' ficld-of-vicw. 

Ground based images of NGC 4631 were also 
recorded in the standard U,B,V, and R band- 
passes and in Ha and H(3 line emission. The ob- 
servations are summarized in Table 1. The FUV, 
Ha, U, B, and l/-band images are shown in Fig- 
ure 1. 

The absolute sensitivity of the UIT was deter- 
mined by observations of stars the spectra of which 
were measured by the lUE (Stecher et al. 1997). 
This procedure resulted in an uncertainty in the 
absolute flux measurements of ~10%. Relative 
sensitivity throughout the Bl passband was mea- 
sured in the laboratory. Relative transmissions 
anywhere in the passband were accurate to 2% 
at transmission levels above 1% of the maximum 
transmission. At 320 nm the instrument response 
is reduced from the maximum response by more 
than 5 dex and at 330 nm by more than 8 dex. 

The UIT image was corrected for geometric im- 
age ( "S" ) distortion according to the formulation 
developed by Greason et al. (1994). Each of the 
ground based images was then transformed to the 
same plate scale and orientation as the UIT image. 

The absolute sensitivity of the U, B, V, and 
R band images was determined using observa- 
tions of the spectrophotometric standard star 
BD-l-33 2642. After establishing the calibration 
using the spectrophotometry of Stone (1977), 
the U, B, and V magnitudes were cross checked 
against the broad band measurements recorded 
by Klemola (1962). The values match within 5%. 

Astromctry of the broad band U,B,V,R and 
H(3 images was determined using stars in common 
with the Digitized Guide Star Catalog of Lasker 
et al. (1990). The astrometric solution of the UIT 
image was based on positions of star forming knots 
in common with the KPNO U-band image yielding 
a positional uncertainty of ± 0.74 arcseconds. 

Continuum subtraction in the Ha image was 
performed using foreground stars in both on-line 



and off-line images. Observations of the planetary 
nebula NGC 7027 were used for absolute calibra- 
tion. Astrometry for the Ha image was adopted 
from Rand et al. (1992), and differs from the 
broad-band astrometry by no more than 1.1 arc- 
seconds. 

The on-line and off-line H(3 images were sep- 
arately calibrated using the spectrophotometric 
standard star, Kopff 27 (Stone 1977). After abso- 
lute calibration, a scale factor of 1.22 was applied 
to the off-line image to allow removal of contin- 
uum sources, i.e. star images from the field. After 
continuum subtraction was performed, the flux of 
the UV-bright region in the disk of NGC 4631, 
CM-67, as measured through a 21" diameter syn- 
thetic aperture was 2.34x10^^*^ ergs s^^ cm^^. 
This value is 6% greater than the equivalent mea- 
surement by Roy et al. (1991). 

3. Morphology 

Inspection of the V and B images shows dust sil- 
houetted against the red stars of the nuclear bulge. 
Most of the dust appears to be concentrated near 
the galactic center; its morphology corresponds 
roughly to the distribution of CO in an interior 
disk as measured by GoUa & Wielebinski (1994), 
Sofue et al. (1990), and of the dust as observed by 
Braine et al. (1995), in 1.3 mm continuum emis- 
sion. Its distribution is chaotic extending above 
the galactic midplane by as much as 1.3 kpc. The 
B, V and Ha images show that the cores of some 
of the CM H II regions are silhouetted against the 
dust patches, implying that much of the bright 
FUV emission originates edgcward from the cen- 
tral 4 kpc of NGC 4631. 

This impression is reinforced by the measured 
H I and H II velocity-position curves. The H I 
velocity-position curves (Weliachew et al. 1978; 
Sofue et al. 1990) imply a constant H I velocity 
at distances greater than ~1' from the galactic 
center. However, at these distances the Ha data 
of Crillon & Monnet (1969) show a linear increase 
in the observed velocity as a function of increas- 
ing distance from the galactic center. Such behav- 
ior would result if the H II gas velocity were con- 
stant beyond ~1' from the galactic center, and 
the emission was from gas located at an approx- 
imately constant galactocentric radius. Since the 
bright FUV/H II regions extend 24 kpc along the 



entire major axis of the imaged galaxy, we deduce 
that most of the H II regions are located near the 
galactic periphery. 

We also note that if the presently observed 
FUV emission regions were to be further obscured 
by dust in an interstellar gas of density 0.5 cm""^, 
as indicated by Rand's (1994) 21 cm measure- 
ments of the NGC 4631 midplanc, an effective gas 
depth of 1.2 kpc would cause more than half of the 
observed UV-bright regions to be undetectable by 
the UIT. To obtain this result, the relationship 

{N{HI+H2)/E{B-V)) = 5.8x10^^ atoms cm-^ mag- 

(Bohlin et al. 1978), determined from measure- 
ments along 75 lines-of-sight in the Galaxy, is as- 
sumed to be valid. We do not include the molecu- 
lar gas density in our estimate since this quantity 
is unknown; its inclusion would further reduce the 
effective gas depth. We also assume that a de- 
tection would be at the 2ct level. Even though 
the dust is observed to be patchy, the fact that 
the entire major axis is delineated by bright FUV 
emission requires that many of the FUV emission 
regions are near the visual periphery of the galaxy. 
Finally we note that peripheral H II regions and 
FUV emitting sources are commonplace in many 
Sc-Sd type galaxies, including MlOl (Waller et 
al. 1997a), M74 (Cornett et al. 1994) and M83 
(Bohlin et al. 1983). 

Generally, the FUV emission coincides closely 
with the Ha emission, much of which shows signif- 
icant departures from the midplane as defined by 
the symmetric distribution of H II regions about 
the major axis (see Crillon & Monnet (1969)). 
If this coincidence is related to the distribution 
of young, massive stars, and contributions to the 
observed FUV light from distant regions of the 
galaxy are masked by the UV opacity of the dis- 
placed gas then star formation must be occurring 
at substantial distances (< 36", 1.3 kpc) above 
the galactic midplane. This does not include the 
even greater extent of Rand's shell 1 perpendicu- 
lar to the galactic midplane (z < 1.7 kpc). Some 
unknown fraction of the high diffuse FUV emis- 
sion at high altitude is likely due to dust scattered 
light. 

Comparison of the FUV and Ha images indi- 
cates that the sizes of the H II and star forming 
regions are about the same as well as being co- 



spatial. This implies that the B-type stars which 
account for most of the FUV light have not moved 
far from their birth place. Typically, a star form- 
ing region which appears to be a single rounded 
area unblended with other regions in the UIT im- 
age subtends an angle of about 8 arcseconds cor- 
responding to a linear distance of 291 pc at the as- 
sumed distance of 7.5 Mpc. Moving with a speed 
of 10 km s~^ stars would disperse to a width of 
291 pc in 14 Myr. Although some of the FUV 
bright regions are clearly blends of smaller regions 
making geometric comparisons impossible, we ex- 
pect that many of the star forming regions must 
be on the order of 14 Myr old or younger. 

4. Photometry 

4.1. FUV and Visual Photometry of 48 
Regions 

Figure 2 is the FUV image which shows hand- 
drawn boundaries and numbering of the apertures 
in which photometry was performed. The criteria 
for establishing the numbered boundaries were (1) 
that each should include at least one unresolved 
FUV emission region that is also seen clearly in 
the B image, (2) that each should include all the 
FUV reflection nebulosity belonging to the se- 
lected emission regions, and (3) that each should 
include all the compact H II regions and Ha fil- 
aments apparently belonging to the selected asso- 
ciations and not to the general background. The 
boundaries thus determined were also applied to 
the U, B, V, Ha and H(3 images. The intent was 
to delineate those data which will yield a consis- 
tent model of star formation within each bound- 
ary. Similarly, a large aperture drawn around the 
entire galactic image in each broad band and emis- 
sion line permitted global quantities to be com- 
puted from spatially integrated fluxes. 

The sky background, subtracted from each im- 
age, is determined from the responses in 50 rect- 
angular apertures distributed over the entire field- 
of-view and located in clear regions between stel- 
lar images. The apertures are 31x31 pixels^ and 
the standard deviation of the distribution of the 
pixel responses within an aperture were taken to 
indicate the rms variation in the pixel-to-pixel re- 
sponse. In flux density units the pixel-to-pixel 
variation is 4.3 x 10^^® ergs cm^^ s^^ A^^ in the 
FUV image and < 1.4 x 10^^^ ergs cm^^ s^^ A^^ 



in the U , B and V images. The sky background is 
the average of the responses in the 50 apertures. 
The standard deviation of the average responses in 
the 50 apertures is taken to indicate the uniformity 
of the sky over the field-of-view. This variation is 
< 5.1x10"^^ ergs cm"^ s"^ A^^ in the FUV im- 
age and < 3.6xl0~^^ ergs cm^^ s^^ A^^ in the C/, 
B and V images. The variation in sky uniformity 
makes a negUgible contribution to the total ran- 
dom error in the derived flux densities. However, 
the contribution of the pixel-to-pixel variation can 
be as large as 28%, although it is generally less 
than 10%. 

Correction for scattered light and sky back- 
ground within the individual apertures was achieved 
in the following manner. Each aperture was ex- 
panded by 2 pixels (2.5 arcseconds) at its edge to 
create annuli between the original and expanded 
apertures. Histograms of the pixel response to the 
background included in the annulus were made 
from which the mean value was taken to be rep- 
resentative of the background from all sources 
underlying the emission from the recently formed 
stars and the gas excited by those stars. The total 
background was simply the mean value multiplied 
by the number of pixels encompassed by the orig- 
inal aperture. 

Other than the correction for the sky back- 
ground, no further correction was made to the 
integrated flux density of each image within the 
large aperture. Expanding and shrinking the large 
aperture by 8 pixels or 20 arcseconds changes the 
integrated flux by 1% or less in all of the broad- 
band images. Thus, the integrated flux is quite 
insensitive to the actual position of the large aper- 
ture on the sky. 

The contribution to the observed FUV contin- 
uum arising from two photon emission in the ion- 
ized gas was evaluated and found to be generally 
at the 1% level with a maximum of 4% of the ob- 
served continuum. This effect was neglected in 
subsequent calculations. 

Error will be made to the extent that the stel- 
lar emission, associated reflected light and line 
emission from the individual regions cannot be 
uniquely isolated from light emitted from the rest 
of the galaxy. This is particularly true for the _B, 
V , and R images and is exacerbated by the edge- 
on aspect of the object. The problem is so severe 
for the _R-band image that the aperture photom- 



etry is very inaccurate, and is for this reason not 
included in the subsequent analysis. In addition, 
uncertainties arising from the background correc- 
tion procedure can be substantial, particularly for 
faint sources. Some indication of the magnitude of 
this systematic error can be had by expanding and 
contracting the apertures by 2 pixels (2.5 arcsec- 
onds). The result is that sky corrected flux den- 
sities change by an average factor of 1.5 for the 
broad-band photometry and factors of 1.4 and 1.7 
respectively for the Ha and Hj3 line photometry. 

Listed in Table 2 are the photometry results 
corrected for background and extinction in the 
Galaxy. The latter correction was accomplished 
using the Galactic value of E(B-V) = 0.007 (de 
Vaucouleurs et al. 1991) (see note b. Table 2) with 
the Seaton (1979) extinction law. The aperture 
number as in Figure 2 and GM H II regions in- 
cluded in each aperture are listed in column 1 and 
2 respectively. The observed FUV (152 nm), [/, B, 
and V fluxes in units of 10~^^ ergs cm~^ s~^ A^^ 
are listed in columns 3, 4, 5 and 6 respec- 
tively, and the Ha and Hj3 fluxes in units of 
10~^^ ergs cm^^ s"^ are listed in columns 7 and 8 
respectively. Values of E{B — V) computed from 
the Balmer Decrement are given in column 9. The 
second number in each column is the la standard 
deviation due to random errors in the measure- 
ments only. The next to last row shows the sums 
of each column for columns 2 through 7, i.e. the 
sum of the fluxes in all 48 apertures. The last row 
shows the flux values integrated over the entire 
galactic image inside the large aperture. 

The B — V color excesses were computed from 
Balmer Decrement and from the difference be- 
tween the measured B — V colors and B — V 
colors determined from models of evolving star 
forming clusters. The models (see Section 5, be- 
low) describe the evolution of clusters in which 
the stars are all coeval, sometimes called "instan- 
taneous burst", or IB models. In principle, the 
derivation of E{B — V) in this way is more real- 
istic than using the Balmer Decrement since it is 
related to the extinction of stellar continuum light 
rather than the extinction of emission from nearby 
excited gas. However, the method suffers from 
the the need for a priori knowledge of the correct 
IB model parameters in order to infer E{B — V). 
With the present data-set, the models were char- 
acterized only by their burst ages which were de- 



termined from the ratio of the number of Lyman 
continuum photons to the FUV luminosity, and 
these latter quantities were determined from the 
FUV and Ha observations corrected for dust ex- 
tinction. Values for all other model parameters 
were assumed as described below in Section 5. 

The distribution oi E{B — V) values determined 
from the Balmer Decrements is shown by the red 
line in Figure 3. A single gaussian fit to the distri- 
bution results in a mean value of E{B — V) = 0.39. 
Though the distribution is sharply peaked, there 
are broad wings at low occurrence levels extend- 
ing to unrealistic negative E{B — V) values. We 
take these broad wings to be the result princi- 
pally of significant errors in the H(3 flux density 
measurements. The less important part of the er- 
ror is statistical; the most important part, as dis- 
cussed above, is due to our inability to identify the 
relevant emission, and is estimated by the vari- 
able aperture photometry. Because of these errors 
we conclude that use of the individual color ex- 
cesses to correct the measured flux densities for 
extinction is inappropriate. Rather, we determine 
a mean value to be applied to all the photome- 
try by means of an iterative procedure. As an 
initial step, we use the mean value of the distri- 
bution computed from the Balmer Decrement to 
correct the FUV and Ha measurements in each 
region for extinction using the extinction law of 
Calzetti et al. (1994). From these corrected values 
we compute models for each region which are char- 
acterized by unique burst times. Subsequently, 
E{B — V) values are derived from the modeled 
and observed colors for each region from which a 
new mean value is determined. This process is re- 
peated until the mean values of E(B^V) converge. 
The final distribution of E{B — V) values, shown 
as a black line in Fig 3, exhibits occurrence lev- 
els in the wings lower than those associated with 
the Balmer Decrement measurements, and a mean 
value of E{B—V) = 0.28 which is less by a factor of 
0.7 than the value derived from the Balmer Decre- 
ment measurements. Qualitatively, this difference 
in derived values of E{B — V) is characteristic of 
stars in regions of active star formation. However, 
the factor by which results of the two methods dif- 
fer is usually smaller, being about 0.44 (Calzetti et 
al. 1994; Fanclh et al. 1988). The 2a width of the 
Gaussian function fitted to the final distribution 
is 0.36 mag. 



Several photometric quantities, which have 
been corrected for internal extinction, are listed in 
Table 3. As in Table 2, the aperture number ap- 
pears in column 1. Magnitudes in the FUV, U, B, 
and V bands are listed in columns 2, 3, 4 and 5 re- 
spectively. The Ha luminosity, the log of the ratio 
of the number of Lyman continuum quanta to the 
FUV luminosity, logio(NLyc/Li52), and the M152- 
U, Mi52-i? and M152-F colors are listed in columns 
6, 7, 8, 9 and 10 respectively. The associated for- 
mal random measurement errors (Icr), excluding 
those arising from the extinction corrections are 
listed adjacent to each quantity in each column. 
Uncertainty in the extinction correction will re- 
sult in significant uncertainty in the listed quan- 
tities. Increasing or decreasing the mean value 
of E{B — V) by the standard deviation in the 
E{B - V) distribution changes the FUV, U, B 
and V magnitudes by 0.80, 0.50, 0.43 and 0.36 
mag respectively, the Ha luminosity by a factor 
of 2 and the M152 - U, M152 - B and M152 - V 
colors by an average of 0.37 mag. The change 
in aperture sizes results in an average change in 
broad-band magnitudes of 0.56 mag, a change in 
the Ha luminosity by 35%, and an average change 
in the M152 - U, M152 - B and M152 - V colors 
of 0.37 mag. In general, random errors in quan- 
tities expressed as magnitudes are comparable to 
the changes induced by variation of the aperture 
sizes. The errors do not include the systematic 
errors associated with the absolute calibration or 
the choice of the reddening law. 

4.2. Integrated Photometry of NGC 4631 

The last line in Table 2 gives the integrated 
photometry of that part of NGC 4631 which can 
be observed, in particular all of the FUV im- 
age inside the large aperture. A comparison with 
two possible analog galaxies, MlOl (Sc(s)I) and 
M83 (SBc(s)II)), both observed face-on is given 
in Table 4. Results from integrated UIT pho- 
tometry of two additional face-on galaxies, M74 
(Sc(s)I) and M51 (Sbc(s)I-II), are also listed. Ex- 
cepting NGC 4631, the FUV luminosities have 
been derived from the results of Waller et al. 
(1997b), with corrections only for extinction aris- 
ing in the Galaxy, and with modifications to take 
into account differences in the galactic distances 
listed by Waller et al. (1997b) and those used 
in this paper. The NGC 4631 152.1 nm flux is 



a factor of 1.1 brighter than the 155.0 nm flux 
(9.81 xlO~^^ ergs cm^^ s^^ A~^) measured by 
Code & Welch (1982) using the University of Wis- 
consin experiment on the first Orbiting Astronom- 
ical Observatory (OAO). The measured NGC 4631 
FUV luminosity is dimmer by factors of 1.2 and 
4.0 than the measured luminosities of M83 and 
MIDI respectively. 

Far infrared (FIR) luminosities after Rice et al. 
(1988), adjusted for the differences in distances be- 
tween that publication and this paper, are listed 
in column 5. The ratios of FIR to FUV lumi- 
nosity are listed in column 6. The comparatively 
high Lfir/^fuv ratio of NGC 4631 compared 
to MlOl, M74 and M51 is expected because of 
its edge-on aspect. If, when viewed face-on, the 
^fir/^fuv ratio for NGC 4631 were similar to 
that of MlOl, i.e '--^4, then we can surmise that 
the FUV luminosity would be -2.3x10^ Lq • If 
the face-on I^fir/^fuv ratio were more nearly 
like that of M74 or M51, then the FUV luminosity 
would be '--^1.1x10^ L0. The exceptional galaxy in 
this group of five is M83 which shows a compara- 
tively large hpjR/hFuv ratio (nearly equal to that 
of NGC 4631) even though it is viewed face-on. It 
may be that the dust in M83 provides especially 
efficient conversion of FUV light into FIR radia- 
tion primarily because of its distribution relative 
to the hot, young stars. The distribution of dust 
in MlOl may also account for the comparatively 
small I^fir/^fuv ratio in that galaxy. 

5. Cluster Ages, Masses, and Star Forma- 
tion Rates 

The photometric results can be combined with 
models of star forming clusters to explore some as- 
pects of the star formation processes in NGC 4631. 
The cluster models we will use are discussed in 
detail in Hill et al. (1994). Here, we will sum- 
marize that discussion. The models are those of 
Landsman (1992), which are based on the work 
of Lequeux et al. (1981) who derive spectral en- 
ergy distributions and masses of evolving star clus- 
ters. The models use the Geneva stellar evo- 
lutionary tracks (Meynet et al. 1994) and the 
models of Kurucz (1991) at logio(Z/Z0) = -0.3. 
In our work, the form of the IMF is ^^(M) — 
dn(M)/d(lnM) ex M"^ for 1.8 Mq < M < 120 M© 
and M-O-e for 0.007 M(d < M < 1.8 M^. We 



assume an upper IMF slope of x = 1.35 conform- 
ing to the Salpeter (1955) IMF. For our purposes, 
an instantaneous burst (IB) model of a cluster 
of co-eval stars is a spectral energy distribution 
(SED), L/b(A, t), corresponding to a given elapsed 
time from a burst of star formation. It is the 
sum of the SEDs of all the stars in the cluster 
at age t weighted by the assumed IMF and nor- 
malized to an astrated mass of 1 Mq. Its units 
are ergs s^^ A~^ Mq^^. The ionizing photon lu- 
minosity is computed by integrating L/b(A, i)/hc 
over A shortward of 912 A. The observed ionizing 
photon luminosity is given by 



N 



Lye 



2.206 LHafi^Ha/hc 



(1) 



where \jHafi is the Ha luminosity obtained from 
the observed flux corrected for both extinction and 
distance. The factor 2.206 is derived under the 
assumption of case B recombination at 10,000K. 
The FUV luminosity, L/s(1520 A,i), is deflned as 
L7b(A, t) integrated over the FUV response curve 
of the UIT. The age of a cluster is simply deter- 
mined by flnding that value of t at which the ob- 
served and model ratios of Nl);c/L/b(1520 A,t) 
are equal. The total astrated mass inferred from 
the IB model is simply the observed luminosity di- 
vided by the normalized IB model SED integrated 
over the response curve for each broad band. 

If continuous star formation (CSF)has occurred 
at an essentially constant rate over a period of 
time, the star formation rate, SFR, can be com- 
puted by integrating the IB models over time; i.e. 
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L, 
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J^LiB{\,t)dt 



MQyr' 



(2) 



where ljobs,\ is the measured spectral luminosity in 
ergs s""'^ A"-'^ and L/b(A, t) is the IB model deflned 
above. The integration in time is over the interval 
of star formation. 

5.1. Instantaneous Burst Models 

Plots of the observed colors, M152-C/, M152- 
B, and Mi52-y, (corrected for extinction), vs. 
logio(N/,j^c/Li52) are shown in Figure 4, a, b, and 
c. Numbers adjacent to the data points in Fig- 
ure 4c identify the aperture corresponding to each 
datum. Error bars indicate the la random photo- 
metric errors. The points which have been plotted 



meet the conditions that the signal to noise ra- 
tios in the M152 — V colors are greater than 3 and 
the color changes due to variations in the aperture 
sizes are less than the la errors. 

The triangle symbols represent the colors com- 
puted from IB models where each triangle corre- 
sponds to a time in the past when the star burst 
occurred in integer units of 1 Myr beginning at 
2 Myrs on the right and ending at 11 Myrs on 
the left. The arrow indicates the direction points 
would move if E(B — V) were to increase by 
0.18 mag which is the la width of the relevant 
E{B — V) distribution (black line) shown in Fig- 
ure 3. All three colors derived from the obser- 
vations of the selected star forming regions tend 
to be more positive than those derived from the 
models. This behavior may be due to our use of 
an inappropriate reddening curve or possibly to 
the assumptions regarding the model input pa- 
rameters. However, in view of the random and 
systematic errors such distinctions are problem- 
atical. For example, increasing the mean value of 
E{B — V) by about 0.11 mag, which is the adopted 
uncertainty in this quantity for integrated values 
as explained below, brings the data in reasonably 
good agreement with the models. A comparison 
of the observational results with the IB model in- 
dicates that if the star formation occurs in bursts 
then the bursts occurred from 2.7 to 10 Myrs in 
the past. It will be noted that the burst times 
are strongly dependent on logio(Nij,c/Li52)- On 
average, the change in aperture size will move the 
points to the left or right an amount equivalent 
to 1 Myr. The implication of the uncertainty in 
the extinction correction is that star formation for 
some regions may have occurred 2.7 Myr further 
in the past, or 2.3 Myr toward the present. 

Individual burst times, determined from the 
measured logio(N/,j,c/Li52) values are computed 
for each aperture. The average variation in burst 
times attributable to the photometric standard de- 
viation is 0.71 Myr, which is about the same as 
variation due to changes in the aperture sizes. The 
total astrated mass is determined for all apertures 
from photometry in the four wide bands as out- 
lined above. The masses summed over the 48 aper- 
tures for each broad band are similar; the largest 
deviation, 20%, is associated with the FUV pho- 
tometry. This result is remarkable in that the 
mass of the stars producing the FUV light consti- 



tutes a small part of the total mass of stars produc- 
ing light in the optical bands, and it implies that 
an average of the summed masses in each band is 
justified. Random errors in the masses summed 
over all the apertures for each band are on the or- 
der of a percent whereas changes in the aperture 
sizes cause the summed masses to vary by a fac- 
tor of 1.41. For integrated quantities corrected for 
extinction internal to NGC 4631 we adopt, some- 
what arbitrarily, an uncertainty in E{B — V) of 
0.11, which is equal to the difference between the 
mean values of the E(B — V) distributions derived 
from the Balmer Decrements and from the stellar 
continua. This difference results in an uncertainty 
in the summed masses by of factor of 2.8. When 
the aperture size variation is combined with the 
extinction uncertainty the summed burst mass is 
uncertain by a factor of 3.9. 

5.2. Continuous Star Formation Models 

Calzetti (1997) has suggested that in active star 
forming regions such as described here, B type 
stars, which supply large amounts of FUV radi- 
ation but little Ha, migrate away from the dust 
clouds of their origin where star formation is con- 
tinuing to produce massive O stars. A typical 
migration time of 50 Myr corresponds to a lin- 
ear dimension of 500 pc and a stellar velocity of 
10 km s~^. Such a picture can explain why the 
color excesses measured from stellar continuum 
colors are often less than those measured from 
Balmer Decrements in the gas surrounding Lyman 
continuum emitting stars. From our observations 
we estimate that the smallest, most compact UV- 
bright regions in NGC 4631 can not have begun 
forming stars substantially more than 14 Myr ago 
purely on the basis of the morphological similarity 
of the Ha and FUV images. The crossing time of 
the largest region is equivalent to a turn-on time 
of 48 Myr in the past. However, the FUV image 
must be composed in part of many spatially un- 
resolved clusters, and with the available data we 
can not disentangle effects due to the diffusion of 
stars from those due to the distribution of stellar 
birth sites. The impression of comparative youth 
is enhanced by the similarity of the mean E(B-V) 
values determine by the Balmer Decrements and 
the integrated stellar continua of each star form- 
ing region. Nevertheless, there is likely to be a 
spread in turn-on times among the UV-bright re- 



gions. The analysis presented in Section 6 implies 
that the large eastern shell may have resulted from 
a star burst which ocurred 20 Myr ago. Somewhat 
arbitrarily, we assume that the most representa- 
tive turn-on time is 14 Myr in the past, but that 
this representative turn-on time could vary by at 
least ± 5 Myr. 

Shown in Figure 4 a, b, and c are colors de- 
termined from CSF models with star formation 
beginning 14 Myrs ago and ending at times indi- 
cated by circles plotted in 1 Myr intervals from 
the present (0 Myrs) on the right to 9 Myrs ago 
on the left. The data points indicate a similar 
range in turn-off times. Star formation turn-off 
times and astration rates computed from relation 
(2) for each passband and the Ha line are listed 
in Table 5, columns 8 through 13. The rate aver- 
aged over all four passbands and the Ha is listed in 
column 14. Omitted from this compilation arc the 
rates for region 28, which shows a turn-off time so 
close to the onset of star formation that the rate is 
unrealistically high. Stated differently, region 28 
must be considered a star burst. The average star 
SFR summed over all the star formation regions is 
0.33 Mq yr^^. If the representative turn-on time 
was 20 Myr ago, then the total SFR would be 
about half this rate, or 0.15 M© yr"-*^. The average 
variation of the summed SFR due to the change in 
the aperture sizes is ±0.17 Mq yr~^, whereas er- 
ror in the extinction correction (0.11 mag) results 
in an uncertainty in the summed SFR of a factor 
of 5.5. The random error in the summed average 
SFR is about 1%. 

The diagnostic diagrams derived for the CSF 
models are very similar to those derived for the 
IB models, their differences being significantly 
smaller than the dispersion of the data. The ob- 
served colors agree with the model colors within 
the limits defined by the estimated errors. The 
slope displayed by the CSF models is flatter and 
less variable than that of the IB models, partic- 
ularly at the most recent turn-off times (or burst 
times). This circumstance arises from the fact that 
most of the FUV flux is provided by B-type stars 
with lifetimes longer than 14 Myrs. Therefore, the 
relative proportion of B and later type stars does 
not vary significantly over the 14 Myr baseline, 
thus yielding approximately constant UV colors 
which vary smoothly with turn- off time. 



6. The Eastern Supershell 

McCray & Kafatos (1987) have developed a 
quantitative picture of the dynamics of super- 
shells, which can be used to interpret our results 
and reinterpret those of Rand & van der Hulst 
(1993). The picture is one in which the supershell 
was formed and driven to its present size primar- 
ily by supernovae occurring in an extraordinarily 
large collection of associations and clusters resem- 
bling those we observe. Combining equations (3) 
and (4) of McCray & Kafatos (1987) gives an ex- 
pression 

i7 = 0.0588(/cm s'^ pc~^) (Rs/Vs) (3) 

where ty, in units of 10 Myr, is the time required 
for the shell to expand to a radius, Rs, with a ter- 
minal velocity of Vs. Adopting Rg = 1.5 kpc and 
Vs = 45 km s~i (Rand & van der Hulst 1993) we 
find tj — 1.96 or about 20 Myr. During this inter- 
val, only those stars with mass equal to or greater 
than M* will become supernovae and contribute 
to the shell expansion. M* is given approximately 
by the relation M* -10 M© (t7/3)-°-^25 (McCray 
& Kafatos 1987), and is 13 Mq for the assumed 
conditions. The total number of stars, N.^, con- 
tributing to shell expansion is given by 

TV* X E51 = no iRs/S2pcf if^ (4) 

where E51 is the energy input to the ISM from 
each supernovae event in units of 10^^ ergs, here 
assumed to be unity, and no is the density of the 
ambient ISM, which is assumed to be uniform. 
Substituting the numerical values for Rg and 17 
gives N* = 2.7x10^ no stars with masses > 13 Mq. 
If no = 0.2 cm^"^, consistent with the estimate of 
Rand & van der Hulst (1993), then N* = 5.3x10'' 
stars (M > 13 Mq). This number is somewhat 
larger than the estimate of Rand & van der Hulst 
(1993) (1 - 3.5) X 10'' stars, and in absolute terms 
is extraordinarily large. For example, Kennicutt 
(1984) has estimated that in 30 Doradus the to- 
tal mass of OB stars between 10 and 100 Mq is 
50,000 Mq. Assuming a Salpeter (1955) IMF, the 
total number of OB stars between 13 Mq and 
100 Mq is 1.54x 10^, a factor of about 34 less than 
our estimated number of OB stars responsible for 
the large shell. We find that the light emitted 
from region 20, which includes the giant H II re- 
gion CM 67, is produced by stars with a total mass 



of 2.7x10'^ Mq. If our assumed IMF is valid, then 
this region contains approximately 1.3x10^ stars 
between 13 M© and 120 Mq, a number which ap- 
proaches that inferred for the ancestral stars of 
the large shell. We infer an even larger number of 
stars (1.5x10"*) for region 28. 

In this picture of shell formation, the UV-bright 
regions associated with the shell, namely those de- 
lineated by apertures 3 through 16, must arise 
from second generation star formation. As Mc- 
Cray & Kafatos (1987) point out, the shell will 
probably burst at high galactic latitudes before 
third generation stars are formed. In fact, as Rand 
& van der Hulst (1993) suggest on the basis of a 
break in the H I emission observed on the north 
side of the shell (their Figure 3), this breakout may 
have already occurred. These authors also point 
out that a shell of 3 kpc diameter can be pro- 
duced initially only if there is significant gas pres- 
sure and/or magnetic pressure at high z. These 
conditions seem plausible when taking into ac- 
count the obvious tidal effects of the interaction 
between NGC 4631 and NGC 4656 manifested in 
the H I emission maps of Weliachew et al. (1978) 
and Rand & van der Hulst (1993). 

If the stars in the vicinity of the supershell were 
made from a reservoir of gas, most of which was 
swept up into the shell, then it is possible to make 
an estimate of a star formation efficiency, here de- 
fined as Msteiiar / (Msteiiar + Mgo^), in the supcr- 
shell region by referring to the Rand & van der 
Hulst (1993) estimate of the hydrogen supershell 
mass. The average total stellar mass in the shell 
region (apertures 3 through 16) inferred from our 
observations in the FUV, U, B, and V bands is 
1.6x10^ M0 based on the IB model. Compari- 
son with the estimated hydrogen shell mass of (1- 
2)xl0^ M© implies an efficiency of -0.047. This 
efficiency represents an upper limit since the con- 
tribution of a molecular gas has not been taken 
into account. 

7. Discussion 

Interpreted in terms of an instantaneous burst 
model, the light in all the numbered apertures 
and in all four wavelength bands implies a total 
created mass of 2.5x10^ Mq within a factor of 
3.9. This number may be underestimated by a 
comparatively small amount due to the absorption 



of light by interstellar grains as discussed below. 
Our estimate can be compared to the equivalent 
number of 5.2x10^ Mq in the Magellanic irregu- 
lar, NGC 4449, which exhibits intense star forma- 
tion (Hill et al. 1994). The difference between the 
two galaxies is even greater than indicated here 
because NGC 4449 is observed nearly face-on re- 
vealing more of its star-forming regions than does 
NGC 4631 which is observed almost edge-on. 

The interior disk of NGC 1068 exhibits nu- 
merous starburst knots surrounding the nucleus. 
Bruhweiler et al. (1991) have estimated the total 
mass of nine of the brightest of these knots, as- 
suming an IB scenario, to be 2.4x10^ Mq, a num- 
ber comparable to that derived for our observed 
48 regions in NGC 4631. Thus, while the masses 
deduced for the star- forming regions of NGC 4631 
are extraordinarily large, they are not uniquely so. 

We also compare the star formation activity 
of the individual star forming regions with the 
prototypical giant H II region, 30 Doradus. Us- 
ing an aperture of 370 pc (d=50 kpc), Kennicutt 
(1984), obtained Lna = 1.5x10^° ergs s"^ and 
an estimated total mass of 50,000 Mq for all stars 
with individual masses between 10 and 100 Mq . 
Extending this Kennicutt estimate over a larger 
mass range by using the Salpeter IMF we find a 
total stellar mass between 0.007 and 120 Mq of 
1.11x10^ Mq. If it is assumed that this mass was 
produced in a burst, we can compare the total 
mass estimate with the IB model mass estimates 
for the individual star forming regions listed in 
columns 3, 4, 5 and 6 of Table 5. Of the 48 
aperture delineated H II regions listed, 6 exhibit 
derived total masses, averaged between the four 
wide-band images, in excess of the 30 Doradus 
value. Taking into account our estimated error, 
this number can be increased to 26. Of these 6 re- 
gions, none exhibit Ha luminosities comparable to 
30 Doradus. Again, the photometric uncertainty 
allows this number to be increased to 3. The Ha 
luminosity of region 20 is about two-thirds that 
of 30 Doradus, but its mass is about a factor of 
2.5 greater. Region 20 incorporates the bright 
H II region CM 67 and is associated with a giant 
molecular cloud, both of which have been tenta- 
tively identified with the tip of a bar (Roy et al. 
1991). Regions 41 and 47 show both Ha luminosi- 
ties and masses comparable to 30 Doradus, and we 
note also that the Ha luminosity of region 41 is 
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comparable to that of the supergiant H II region 
in MlOl (NGC 5461), which has an Ha luminos- 
ity uncorrected for extinction of 2.7x 10^° ergs s~^ 
(Waller 1990). However, with some significant ex- 
ceptions such as those highlighted above, most of 
the NGC 4631 star forming regions are less mas- 
sive and older than 30 Doradus. Even when the 
upper limit on E{B — V) is used to correct for in- 
ternal extinction only 7 regions exhibit Ha lumi- 
nosities about equal to or exceeding that of 30 Do- 
radus. 

A "short term" star formation rate (SFR) av- 
eraged over the range in burst times exhibited by 
the data, i.e. 7.3 Myrs, can be computed using the 
created masses summed over all 48 apertures and 
averaged over the four broad bands. This rate is 
~3.5 Mq yr~^ which is substantial. As discussed 
below, it is about equal to the global SFR deter- 
mined from FIR measurements, and a factor of 7 
less than the total SFR measured by Dahlem et 
al. (1995). 

In summary, when the FUV and optical im- 
agery is interpreted in terms of the instantaneous 
burst paradigm it is obvious that NGC 4631 incor- 
porates exceptionally large star formation activity. 

Interpretation of our results in terms of CSF 
models points to a galaxy of a more ordinary va- 
riety. The total average star formation rate deter- 
mined from the FUV data summed over the 48 
apertures for NGC 4631 is 0.33 Mq yr^^, which 
is a factor of 2.7 less than than the equivalent 
NGC 4449 rate (0.88 M© yr"!). However, our as- 
sumed uncertainty in the mean value of E{B — V) 
allows our derived SFR to exceed the NGC 4449 
rate by a factor of 2. Thus, that part of NGC 4631 
which is observable and has incorporated star for- 
mation within the past 14 Myr creates stars at 
a rate comparable to the total rate of the face- 
on starburst galaxy, NGC 4449. Assuming that 
the turn-on time was 20 Myr ago, which results in 
an SFR equal to 0.15 M© yr~^, does not change 
this conclusion. The NGC 4631 results can also be 
compared to the high rate of 2.4 M© yr~^ found by 
Bruhweiler et al. (1991) when interpreting in terms 
of a CSF scenario their observations of the star 
forming regions found in the disk of NGC 1068. 
Assuming constant star formation, our observa- 
tions indicate a total mass created in 14 Myr of 
about 4.6x10^ Mq, which is less than the mass 
inferred from the IB model (5.1x10^ M©) by a 



factor of 5.5. This factor increases to 8.4 if the 
turn-on time was 20 Myr ago. 

Dahlem et al. (1995) have related radio contin- 
uum observations of the halos of several edge-on 
galaxies, including NGC 4631, to the rates of su- 
pernovae occurring within the galactic disks. The 
inferred supernova rate for NGC 4631 is 0.23 yr^^. 
Assuming equilibrium between the birth rate and 
death rates of supernovae producing stars and a 
Salpeter (1955) IMF over the mass range between 
0.007 and 120 M© we find that this result implies 
a star formation rate of 24 M© yr~^. 

A significantly smaller rate is derived from the 
FIR luminosity. For stars with masses greater 
than 10 M0 



SFRpmiMQ yr-^) - 1.4 x lO-^^LpmiLQ) (5) 

(Sauvage & Thuan 1992). In the case of NGC 4631, 
SFRi./K(M0 > 10) = 1.3 Mq yr-i. Employing 
the Salpeter IMF described in Section 5, the SFR 
for stellar masses between 0.007 Mq and 120 Mq 
is 8.1 Mq yr-i or 4.5 M© yr'^ if 45% of the FIR 
luminosity is from the diffuse interstellar medium. 

The total continuous SFR, summed over the 
48 apertures in our data set, 0.33 M© yr~^, is far 
smaller than the rate inferred from Dahlem's radio 
observations. It, of course, refers only to the star 
forming activity delineated by the 48 apertures. 
If it is assumed that the recorded FUV emission 
from the whole galaxy arises from star formation, 
then multiplication of the 0.33 Mq yr~^ rate by 
the ratio of the observed FUV flux in the large 
aperture to the observed flux summed over all 48 
small apertures increases the continuous SFR to 
2.2 Mq yr^^, still a very small SFR compared to 
the Dahlem value but one that approaches the FIR 
rate. 

It is possible that a large portion of the newly 
created stars are so heavily shielded by dust that 
there is no indication of their existence in any 
of our wavelength bands. For example, Bruh- 
weiler et al. (1991) estimate that the emerging 
UV light from FUV bright knots in the disk of 
NGC 1068 implies an SFR which is only 7% of 
the true value. This estimation was based on the 
relative FIR and UV knot luminosities and the as- 
sumption that most of the FIR emission is repro- 
cessed UV radiation from hot stars absorbed and 
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reradiated by interstellar grains. Unfortunately, 
no far infrared observations of NGC 4631 exist 
at a spatial scale which would allow the isolation 
of FIR emission from the galactic background for 
any of the star formation regions. An alternative 
procedure for approximating the relative contri- 
bution of reprocessed UV radiation in NGC 4631 
star-forming regions is to measure the ratio of in- 
tegrated FUV and FIR luminosities of face-on 
galaxies. The galaxy M74 (NGC 628) is ideal for 
this purpose since the required luminosities are 
known and E{B — V) =0.33 over much of the 
galaxy (Cornett et al. 1994), a value which is sim- 
ilar to the mean value of the color excesses mea- 
sured for the NGC 4631 star-forming regions. The 
ratio of the total integrated FIR luminosity from 
IRAS measurements (Rice et al. 1988) to the inte- 
grated FUV luminosity from UIT measurements 
(Waller et al. 1997b) is 8.5 (see Table 4). Similar 
ratios for other face-on galaxies are also listed in 
Table 4. 

It is likely that the dust emitting the IR radia- 
tion is heated in part by the interstellar radiation 
field, with an approximate contribution to the to- 
tal IR heating budget of 45% (Sauvage & Thuan 
1992); though see contrary arguments of Devereux 
& Young (1992) and Devereux & Scown (1994) 
who claim that almost all of the IR radiation is 
reprocessed FUV emission. Assuming that such 
processes are occurring in NGC 4631, and adopt- 
ing the M74 FIR/FUV luminosity ratio of 8.5 
(the average FIR/FUV luminosity ratio for the 
face-on spiral galaxies M83, MlOl, M74, and M51 
is 8.8) and an interstellar to total IR radiation 
energy density ratio of 0.45, we find that the SFR 
determined from the extinction corrected FUV lu- 
minosity will be increased by a factor of only 1.21 
when taking into account the stellar FUV light ab- 
sorbed by the interstellar grains. Thus, the total 
SFR for the 48 star forming regions alone becomes 
0.40 Mq yr""'^ with a large uncertainty. If the inte- 
grated FUV emission within the large aperture is 
assumed to arise entirely from star formation then 
the total estimated SFR in that part of NGC 4631 
which we are able to observe becomes 2.7 Mq yr"-'^. 

The correction is not sensitive to assumptions 
about the IR heating budget. For example, if all 
the observed FIR light was reprocessed FUV light 
from star formation instead of only 55% as we have 
assumed then the correction factor would increase 



by only 16%. On the basis of observations of M83 
and MlOl the uncertainty in the value of R, i.e. 
the applicable value of the ratio of observed FIR 
to FUV light, induces an uncertainty in the cor- 
rection factor of 12%. However, the correction 
is still quite uncertain since FIR observations of 
NGC 4631 with spatial resolution comparable to 
those of the UIT have not been made. 

If all the FIR light is reprocessed FUV light, 
then the SFR inferred from the FUV light in the 
large aperture becomes 3.1 Mq yr~^, again with 
large uncertainties, which should be compared to 
the SFR computed from the observed FIR emis- 
sion via relation (5) of 8.1 Mq yr^^. The result 
that the SFR summed over the 48 apertures is 
from factors of 11 to 25 less than the SFR deter- 
mined from the FIR observations is not surprising 
considering the edge-on aspect of NGC 4631. Al- 
though the amount of created stellar mass differs 
by a factor of 5.5 there is no reason for distinguish- 
ing between the applicability of the IB and CSF 
models on the basis of the diagnostic diagrams in 
Figure 4. The fact that the large aperture FUV 
photometry leads to an SFR about equal to that 
derived from the FIR within a factor of 2 is sur- 
prising, and may indicate that most of the FUV 
light arises from short term star bursting activ- 
ity, i.e. that the exclusive application of the CSF 
model is not appropriate. 

By contrast, the SFR inferred from the radio 
halo (Dahlem et al. 1995) is a factor of at least 5 
larger than that associated with either the large 
aperture FUV photometry or the FIR results (as- 
suming only 55% of the observed FIR is repro- 
cessed FUV light). We speculate that the SFR 
derived from the radio observations is too large be- 
cause the assumption of equilibrium between the 
birth and death rates of the supernovae producing 
stars is invalid, and that the radio halo is due to 
a starburst and not continuous star formation. 

Finally, to the extent that the IB paradigm is 
valid, our conclusions regarding the formation of 
the supershell in NGC 4631 indicates that strong 
starburst activity has occurred in that part of 
NGC 4631 approximately 20 Myrs ago. 

8. Summary 

A far ultraviolet image of the galaxy NGC 4631 
centered on 152.1 nm was recorded by the Ultra- 
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violet Imaging Telescope. The FUV image, com- 
bined with similar images in the U, B, V, and R 
bands and in Ha and H(3, is used to study star for- 
mation in this late type (Sc/SBd) galaxy, viewed 
nearly edge-on, with the following results. 

1. For a disk galaxy seen edge-on, NGC 4631 is 
extraordinarily bright at FUV wavelengths. The 
attenuation of light from 48 star forming regions 
located uniformly along the whole length of the 
visible galactic disk is low with the mean value of 
E{B — V) =0.28. Dust is seen silhouetted against 
the cool star population, and the star forming 
regions appearing in Ha emission are seen sil- 
houetted against the dust. These facts together 
with our interpretation of previously measured H I 
and H II position-velocity curves indicate that the 
FUV emission arises near the edge of the observ- 
able galaxy, as is observed in the face-on Sc/SBc 
galaxies MlOl and M83. 

The overall FUV morphology of this edge-on 
disk galaxy is remarkably similar to that of the 
so-called chain galaxies that have been recently 
observed at high redshift, thus implying a similar 
interpretation for at least some of those distant 
FUV sources (see Smith ct al. (1997)). 

2. Photometric data is presented in the form 
of diagnostic diagrams showing the M152 — U, 
M152 — B and M152 — V colors plotted against 
logio(Nz,yc/Li52)- The data exhibit considerable 
dispersion due to errors made in the determination 
of the relevant background and extinction correc- 
tions. Both instantaneous burst (IB) models with 
burst ages ranging from 2.7 to 10 Myr and continu- 
ous star formation (CSF) models with star forma- 
tion beginning 14 Myr ago and terminating from 
to 9 Myr ago fit the data equally well. The uncer- 
tainties permit the burst ages and turn-off times 
to be extended 1 to 3 Myrs to earlier times. 

3. The ratio of the mean value of E{B — V) de- 
termined from the continua of stars to that found 
from measurements of Balmer Decrements is 0.7, 
and is thus greater than the value of 0.44 quoted 
by Calzetti (1997) as being typical of regions ex- 
hibiting active star formation. This result is com- 
patible with our conclusion on the basis of FUV 
emission morphology that the observed star for- 
mation is young, i.e. that the stars have not had 
much time to migrate from their birth site. 

4. If the observational results are interpreted 



with the IB models, the inferred stellar masses 
created in each star forming region are substan- 
tial. As determined by our analysis, the inferred 
stellar masses of 6 regions exceed the mass of 
30 Doradus. Errors permit this number to in- 
crease or decrease by up to 20. The region of 
the second largest derived mass, 2.7x10^ Mq , 
incorporates the bright H II region CM 67 and 
is associated with a giant molecular cloud, both 
of which have been tentatively identified with the 
tip of a bar (Roy et al. 1991). Of these 6 re- 
gions, none exhibit Ha luminosities equal to that 
of 30 Doradus. Photometric uncertainties permit 
this number to increase to 3. The Ha luminosity 
of region 41 is comparable to that of the super- 
giant H II region in MlOl (NGC 5461). The total 
created stellar mass deduced from the light emit- 
ted from the star forming regions is 2.6x10^ Mq 
to within a factor of 2.3, which is about equal to 
the stellar mass (2.4x10^ Mq) of the UV emitting 
knots in the disk of NGC 1068 (Bruhweiler et al. 
1991) and is a factor of 5 greater than the stel- 
lar mass (5.2x10^ Mq) associated with the FUV 
emitting knots in NGC 4449 (HiU et al. 1994). Us- 
ing the range in burst times of 7.3 Myrs, we de- 
rive an average "short term" star formation rate 
of 3.5Mq yr^i. 

5. Based on previously obtained data, appli- 
cation of the shell generation theory of McCray 
& Kafatos (1987) shows that the number of su- 
pernova producing stars which could generate the 
3 kpc diameter ring of FUV and H I emission 
near the east end of the galactic image is 5.3 x lO"', 
which is somewhat greater than the range (1 to 
3.5) X 10** estimated by Rand & van der Hulst 
(1993) based on the dynamics of a hydrogen su- 
pershell surrounding the stellar ring. While this 
number is extraordinarily large, it is comparable 
to the same number, 1.3x10^, associated with the 
intense star forming region delineated by aper- 
ture 20, a region that includes the giant H II re- 
gion CM67. Those FUV bright regions on and 
within the periphery of the large shell (apertures 
3 through 16) must represent second generation 
star formation, and this conclusion together with 
the shell analysis implies that strong star forma- 
tion occurred approximately 20 Myrs ago in that 
part of NGC 4631 incorporating the large shell. 

6. Analysis of the FUV, U, B, and V photome- 
try in the context of the CSF paradigm yields a to- 
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tal SFR for the 48 individual star forming regions 
equal to 0.33 M0 yr~^, uncertain by a factor of 
8.2. Approximate accounting for UV light which is 
absorbed by interstellar grains leads to a corrected 
rate of 0.40 Mq yr"-'^. Under the assumption that 
55% of the observed FIR is from reprocessed FUV 
light this corrected rate is about a factor of 11 
less than the SFR inferred from the measured to- 
tal FIR (4.5 M0 yr~^), which is not surprising 
considering the edge-on aspect of NGC 4631. If 
the observed FUV light integrated over that part 
of the galaxy we are able to observe is assumed 
to arise in constant star formation, then the SFR 
equals 2.6 M© yr"-'^ when approximate accounting 
for the reprocessing of FUV light by interstellar 
grains is taken into account. This is a surprising 
result in that one would expect the SFR inferred 
from the integrated FUV data to be significantly 
less than the global SFR determined from the FIR 
observations, and the photometric uncertainties 
exacerbate the problem. We interpret this para- 
dox to indicate that the FUV emission should not 
be seen exclusively in terms of the CSF model, but 
that much of the observed global FUV light arises 
in star bursts. This interpretation is consistent 
with the inference of star bursting activity near 
the galactic center drawn by GoUa & Wielebinski 
(1994) and in an inner disc (r«50") by Duric et 
al. (1982). 
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Fig. 1.— Far-ultraviolet (FUV), Ha, U, B, and V, images of NGC 4631. North is up, east to the left. 
Portions of the dwarf elliptical galaxy, NGC 4627, which is part of the interacting system composed of 
NGC 4631, NGC 4656 and NGC 4627, can be seen on the north side of the U, B, and V images. 
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Fig. 2. — FUV image of NGC 4631 as in Figure 1 with aperture boundaries indicated and numbered. 
The apertures are numbered, 1 to 48, from left to right, but because of crowding only the units number 
is shown except for apertures numbered 10, 20, 30 and 40. The large aperture boundary encompassing 
the whole image is used to compute global properties. 
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Fig. 3. — Histograms showing distribution of E{B — V) values of each star forming region. Red line: 
distribution based on the measured Balmer Decrements. Black line: distribution based on the difference 
between modeled colors and observed colors. Single Gaussian fits yield mean values of E{B — V) equal 
to 0.39 mag and 0.28 mag respectively, indicated by vertical dashed lines. The 2ct width of the Gaussian 
function fitted to the distribution represented by the black line is 0.36 mag. 
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Fig. 4a. — Color vs. logio(Nij,c/Li52) diagrams. Colors M152-U, M152-B and M152-V arc plotted in 
panels a, b and c respectively. Error bars indicate la errors. Numbers in panel c identify associated 
apertures. Filled circles: observational results corrected for dust attenuation using the formalism of 
Calzetti et al. (1994). Open triangles: values derived from instantaneous burst models computed at 
integer burst times beginning at 2 Myr on right and increasing to 11 Myr on left. Open circles: values 
derived from continuous star formation models computed at integer turn-off times beginning at Myr 
on right and increasing to 9 Myr on left. The arrow indicates the effect of increasing the correction for 
extinction on the observed colors arising from an increase in E{B — V) of 0.18 mag. 
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Fig. 4b.— Mi52-B vs. login (Nlj,c/Li52) diagram. 



21 



-5 



> 




12.0 



12.5 

log(NLy 



13.0 



13.5 



c/ ^152; 



Fig. 4c.— Mi52-y vs. login (NLyc/Li52) diagram. 
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Table 1 
Observations 



Instrument 


Date 


Exposure 
(seconds) 


Bandpass 
Ao/AA 


Resolution 

(FWHM) 


UIT 


03/10/1995 


1140 


FUV, 152.1/35.4 nm 


3.4" 


Kitt Peak 0.9m 


05/16/1996 
05/16/1996 
05/16/1996 
05/16/1996 


600 
300 
300 

180 


U 
B 
V 
R 


3.0" 
1.6" 
1.1" 
1.6" 


Palomar'^1.5m 


05/1990 


2000 


Ua, 658.5/1.5 nm 


1.2" 


Mt. Laguna 


05/23/1996 


600 


H/3, 487.8/5.0 nm 


4.5" 



^The Ha data was provided by Dr. R. J. Rand. See Rand et al. (1992). 



Table 4 
FUV AND IR Luminosities of Galaxies 



Name 


AB^ 


Distance^ 


Lfc/v'' 


LlR^ 


^ir/^fuv 




(Mag) 


(Mpc) 


Lq 


Lo 




NGC 4631 


0.02 


7.5 


6.99x10^ 


9.34x10^ 


13.40 


NGC 5236 (M83) 


0.14 


4.8 


8.68x10^ 


1.20x10^° 


13.80 


NGC 5457 (MlOl) 


0.12 


7.7 


2.77x10^^ 


1.12x10^° 


4.04 


NGC 628 (M74) 


0.12 


7.8 


3.34x10^ 


2.85x10^ 


8.53 


NGC 5194 (M51) 


0.00 


8.4 


1.23x10^ 


1.08x101° 


8.78 



^Foreground Galactic extinction: NGC 4631, NGC 628, NGC 5194, Burstein 
and Heiles 1984; NGC 5236, AND NGC 5457, Hill et al.(1994). 

•^NGC 5236, Pierce (1993); NGC 5457 and NGC 5194, Feldmeir and CiarduUo 
(1997); NGC 628, Sharina (1996). 

'^Except for NGC 4631, luminosities are computed from Waller (1997b), and 
are adjusted for the distances given in Table 4. 

'ilRAS luminosities are from (Rice et al. 1988) and are adjusted for distances 
given in Table 4. 
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TABLE 2 

Observations Corrected for Background and Galactic Extinction 



Ap. 


CM. 


p ab 
1152 


f„-' 




fs"" 


fv"' 




fifo 


be 


tH0 


,c 


E(B-V)' 


1 




ID 


(flux) 


('^) 


(flux) 


i'^) 


(flux) 


i<y) 


(flux) 


(<^) 


(flux) 


M 


(flux) 


M 


(excess) 


M 


(1) 


(2) 


(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




1 


88 


1.50 


0.15 


1.42 


0.05 


1.47 


0.03 


1.03 


0.02 


193.39 


1.80 


109.85 


3.25 


-0.44 


0.03 


2 


87 


1.84 


0.14 


0.86 


0.05 


0.81 


0.04 


0.60 


0.03 


108.01 


1.27 


38.37 


4.28 


-0.02 


0.10 


3 


86 


0.95 


0.10 


0.25 


0.03 


0.08 


0.03 


0.05 


0.02 


50.75 


1.28 


21.56 


3.04 


-0.18 


0.13 


4 


85,83 


7.26 


0.26 


5.07 


0.11 


4.31 


0.06 


3.06 


0.07 


272.19 


4.38 


27.22 


6.64 


1.14 


0.22 


5 


84 


3.66 


0.17 


2.30 


0.05 


2.32 


0.05 


1.73 


0.05 


114.96 


2.02 


26.52 


4.89 


0.38 


0.17 


6 




4.17 


0.21 


2.05 


0.08 


2.03 


0.09 


1.34 


0.06 


22.58 


2.24 


8.87 


4.02 


-0.11 


0.42 


7 


81 


7.15 


0.22 


4.05 


0.06 


3.06 


0.06 


1.98 


0.04 


666.94 


3.92 


201.98 


5.19 


0.13 


0.02 


8 


79 


0.39 


0.14 


0.21 


0.04 


0.34 


0.04 


0.18 


0.03 


23.54 


1.13 


15.97 


3.62 


-0.61 


0.21 


9 


80 


1.05 


0.11 


0.61 


0.04 


0.74 


0.03 


0.49 


0.04 


148.76 


2.71 


45.50 


3.40 


0.12 


0.07 


10 


77b 


3.57 


0.17 


1.37 


0.07 


1.10 


0.07 


0.69 


0.05 


97.70 


2.36 


13.50 


3.07 


0.84 


0.21 


11 


77 


0.86 


0.11 


0.48 


0.04 


0.31 


0.04 


0.19 


0.03 


49.79 


3.08 


7.55 


3.60 


0.76 


0.44 


12 


78 


2.84 


0.20 


1.29 


0.05 


1.14 


0.06 


0.82 


0.05 


61.31 


1.32 


13.27 


3.34 


0.44 


0.23 


13 




0.50 


0.12 


0.15 


0.04 


0.49 


0.02 


0.16 


0.02 


29.69 


1.16 


3.41 


2.74 


1.01 


0.73 


14 


71,74 


9.72 


0.29 


2.93 


0.11 


2.70 


0.13 


1.65 


0.09 


81.83 


1.90 


106.67 


6.13 


-1.20 


0.06 


15 


73 


0.67 


0.09 


0.49 


0.04 


0.43 


0.03 


0.49 


0.03 


49.20 


0.88 


18.59 


3.21 


-0.07 


0.16 


16 


72 


0.17 


0.10 


0.21 


0.03 


0.18 


0.03 


0.21 


0.02 


75.12 


1.52 


9.69 


3.31 


0.91 


0.31 


17 




0.68 


0.16 


0.90 


0.06 


0.81 


0.06 


0.72 


0.05 


9.46 


1.23 


7.85 


4.56 


-0.79 


0.54 


18 


70 


1.43 


0.17 


0.70 


0.06 


0.82 


0.07 


0.58 


0.05 


37.35 


0.83 


7.50 


3.62 


0.50 


0.44 


19 


68,69 


2.04 


0.27 


1.41 


0.11 


1.64 


0.12 


0.60 


0.09 


40.08 


2.68 


31.00 


3.69 


-0.73 


0.12 


20 


65,66,67 


13.46 


0.45 


9.15 


0.24 


8.71 


0.28 


6.90 


0.21 


524.85 


7.52 


118.44 


5.78 


0.40 


0.05 


21 


64 


1.75 


0.18 


0.85 


0.06 


0.91 


0.07 


0.78 


0.06 


51.98 


4.09 


13.01 


3.42 


0.30 


0.25 


22 


59,61,63 


7.50 


0.25 


2.66 


0.10 


0.93 


0.27 


0.54 


0.31 


323.15 


4.86 


93.46 


5.84 


0.17 


0.06 


23 


54,55,58 


1.89 


0.18 


1.40 


0.08 


1.44 


0.12 


1.82 


0.11 


168.61 


3.15 


37.68 


4.44 


0.41 


0.11 


24 


56,57 


2.11 


0.17 


1.44 


0.09 


2.32 


0.12 


1.87 


0.10 


167.23 


3.66 


15.04 


4.64 


1.24 


0.28 


25 


52 


1.17 


0.19 


1.00 


0.06 


1.19 


0.03 


0.63 


0.07 


125.63 


3.78 


5.29 


3.61 


1.93 


0.62 


26 


47,48 


2.93 


0.26 


1.01 


0.07 


0.90 


0.07 


0.55 


0.05 


117.71 


3.89 


37.32 


3.12 


0.09 


0.08 


27 


45,46 


5.07 


0.17 


1.70 


0.06 


0.81 


0.06 


0.44 


0.05 


348.28 


4.64 


63.78 


4.10 


0.59 


0.06 


28 




2.72 


0.17 


2.41 


0.08 


2.77 


0.09 


2.44 


0.07 


5.95 


2.46 


68.83 


3.04 


-3.19 


0.38 


29 


43 


2.43 


0.17 


1.07 


0.07 


1.65 


0.08 


0.94 


0.06 


236.39 


4.87 


58.61 


4.16 


0.31 


0.07 


30 


42 


2.53 


0.18 


1.36 


0.08 


1.24 


0.09 


0.90 


0.07 


34.74 


1.52 


10.70 


2.98 


0.11 


0.26 



:0003 das 6 I Age 1 6000/^^-0 jjs^:atxj^: 



TABLE 2— Continued 



Ap. 


CM. 


fl52-" 


fu^' 


= 


ffl" 


b 


fv" 


b 


fH.'"' 




fH,'^ 


c 


E(B-V) 


d 




ID 


(flux) 


i'^) 


(flux) 


M 


(flux) 


M 


(flux) 


M 


(flux) 


M 


(flux) 


i'^) 


(excess) 


i-y) 


(1) 


(2) 


(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




31 


41 


2.07 


0.20 


0.24 


0.05 


0.17 


0.05 


0.02 


0.03 


130.26 


5.47 


35.14 


3.56 


0.24 


0.10 


32 


40 


3.89 


0.30 


1.81 


0.09 


1.75 


0.09 


0.96 


0.07 


94.63 


3.63 


7.11 


4.21 


1.40 


0.54 


33 


39 


0.94 


0.12 


0.30 


0.03 


0.25 


0.03 


0.14 


0.02 


63.34 


2.67 


14.96 


2.61 


0.36 


0.16 


34 


37,38 


14.51 


0.47 


5.24 


0.12 


4.42 


0.12 


2.70 


0.08 


562.53 


8.19 


109.66 


4.89 


0.53 


0.04 


35 


35,36 


2.81 


0.27 


1.27 


0.07 


1.31 


0.07 


0.86 


0.05 


107.94 


2.41 


4.28 


5.56 


1.98 


1.19 


36 


34 


6.64 


0.36 


2.52 


0.12 


1.83 


0.14 


1.22 


0.09 


487.84 


10.70 


87.00 


4.20 


0.61 


0.05 


37 


31,33 


2.69 


0.23 


1.23 


0.08 


0.83 


0.06 


0.50 


0.04 


309.57 


14.31 


67.26 


5.25 


0.43 


0.08 


38 


27,29 


5.55 


0.24 


1.75 


0.07 


1.58 


0.07 


0.67 


0.05 


181.37 


6.09 


27.01 


3.61 


0.78 


0.13 


39 


28 


1.96 


0.18 


0.91 


0.06 


0.95 


0.05 


0.72 


0.03 


182.95 


6.44 


21.64 


3.80 


0.99 


0.16 


40 


26 


1.05 


0.12 


0.00 


0.04 


0.38 


0.03 


0.20 


0.03 


28.81 


1.83 


16.19 


3.54 


-0.43 


0.21 


41 


23,24 


12.10 


0.26 


5.20 


0.08 


3.53 


0.06 


2.47 


0.04 


1023.30 


11.98 


223.14 


5.69 


0.43 


0.03 


42 


22 


1.40 


0.14 


1.01 


0.05 


0.71 


0.04 


0.49 


0.03 


256.29 


6.05 


61.54 


3.58 


0.34 


0.06 


43 


20,21 


0.43 


0.14 


0.49 


0.05 


0.35 


0.04 


0.21 


0.02 


76.73 


3.20 


11.75 


3.03 


0.75 


0.24 


44 


19 


4.51 


0.21 


2.40 


0.06 


1.98 


0.05 


1.27 


0.03 


368.20 


8.94 


76.72 


3.54 


0.47 


0.05 


45 


15,16 


0.66 


0.12 


0.27 


0.04 


0.27 


0.03 


0.18 


0.02 


80.43 


7.34 


4.46 


2.95 


1.68 


0.61 


46 


14 


1.54 


0.12 


0.95 


0.04 


0.63 


0.02 


0.49 


0.02 


192.46 


5.94 


26.80 


3.03 


0.84 


0.11 


47 


7,8,9 


12.51 


0.30 


5.75 


0.10 


3.79 


0.09 


2.41 


0.06 


730.31 


8.27 


165.52 


6.14 


0.39 


0.04 


48 


1,3 


1.47 


0.14 


0.46 


0.04 


0.47 


0.03 


0.18 


0.02 


64.47 


2.11 


18.03 


3.92 


0.20 


0.20 


Sum^ 




170.71 


1.48 


82.58 


0.54 


72.83 


0.61 


50.04 


0.52 


9178.61 


35.40 


2215.25 


28.84 


0.00 


0.00 


N4631' 




1121.83 


3.71 


677.66 


0.93 


790.56 


0.51 


680.17 


0.37 


25328.77 


3.97 


823.65 


36.58 


0.00 


0.00 



^Units: 10"^^ ergs sec ^ cm ^ A ^ 

Used Galactic extinction curve from Scaton (1979) and A^ — -03 from dc Vaucoulcurs ct al. (1991) to derive E(B-V) — A^ / 4.1 — 0.007 
'^Units: 10"" ergs see" cm~ 

E(B-V) computed from Balmcr decrement 
■^Sum of non-negative fluxes in the apertures 

Total flux in large aperture 



TABLE 3 

Observed Magnitudes and Colors Corrected for Internal Extinction 



Ap. 


152' 




U 




B 




V 




Life 


b,c 


log(Ni,„e / 


Ll52)'^ 


mi52 - 


mu 


mi52 - 


m_B 


mi52 - 1 


tnv 




(mag.) 


M 


(mag.) 


M 


(mag.) 


M 


(mag.) 


i'^) 


(lum.) 


i'^) 


(ratio) 


i'^) 


(mag) 


i'y) 


(mag) 


M 


(mag) 


{•y) 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




(10) 




1 


12.86 


0.38 


14.27 


0.23 


14.97 


0.19 


14.99 


0.16 


3.85 


0.49 


13.20 


0.05 


-1.41 


0.17 


-2.11 


0.19 


-2.13 


0.22 


2 


12.63 


0.93 


14.81 


0.58 


15.62 


0.50 


15.57 


0.42 


2.15 


0.70 


12.86 


0.09 


-2.17 


0.34 


-2.99 


0.41 


-2.93 


0.49 


3 


13.35 


7.20 


16.17 


4.48 


18.19 


3.88 


18.37 


3.24 


1.01 


2.54 


12.82 


0.64 


-2.82 


2.45 


-4.84 


3.02 


-5.02 


3.62 


4 


11.14 


0.28 


12.89 


0.18 


13.80 


0.15 


13.81 


0.13 


5.41 


0.54 


12.66 


0.04 


-1.74 


0.13 


-2.65 


0.16 


-2.66 


0.19 


5 


11.89 


0.36 


13.75 


0.22 


14.47 


0.19 


14.43 


0.16 


2.29 


0.29 


12.59 


0.05 


-1.86 


0.17 


-2.58 


0.20 


-2.54 


0.24 


6 


11.75 


0.89 


13.87 


0.56 


14.61 


0.48 


14.70 


0.40 


0.45 


0.15 


11.82 


0.09 


-2.12 


0.30 


-2.87 


0.37 


-2.95 


0.44 


7 


11.16 


0.41 


13.13 


0.26 


14.17 


0.22 


14.28 


0.18 


13.26 


1.90 


13.06 


0.04 


-1.97 


0.14 


-3.01 


0.17 


-3.12 


0.20 


8 


14.32 


2.98 


16.35 


1.85 


16.55 


1.60 


16.87 


1.33 


0.47 


0.48 


12.87 


0.28 


-2.03 


0.98 


-2.23 


1.14 


-2.55 


1.35 


9 


13.25 


0.89 


15.19 


0.55 


15.71 


0.47 


15.80 


0.40 


2.96 


0.91 


13.24 


0.11 


-1.94 


0.42 


-2.47 


0.50 


-2.55 


0.60 


10 


11.92 


1.25 


14.31 


0.78 


15.28 


0.67 


15.43 


0.56 


1.94 


0.85 


12.53 


0.11 


-2.40 


0.41 


-3.37 


0.50 


-3.51 


0.60 


11 


13.47 


2.51 


15.44 


1.56 


16.64 


1.35 


16.83 


1.13 


0.99 


0.87 


12.85 


0.23 


-1.97 


0.85 


-3.17 


1.04 


-3.36 


1.24 


12 


12.17 


0.95 


14.37 


0.59 


15.24 


0.51 


15.23 


0.43 


1.22 


0.40 


12.42 


0.10 


-2.21 


0.37 


-3.08 


0.45 


-3.06 


0.54 


13 


14.05 


1.68 


16.69 


1.07 


16.16 


0.89 


17.04 


0.75 


0.59 


0.34 


12.86 


0.21 


-2.64 


0.77 


-2.11 


0.86 


-2.98 


1.03 


14 


10.83 


0.99 


13.48 


0.61 


14.31 


0.53 


14.48 


0.44 


1.63 


0.56 


12.01 


0.09 


-2.65 


0.33 


-3.48 


0.40 


-3.65 


0.48 


15 


13.73 


1.00 


15.42 


0.62 


16.31 


0.54 


15.79 


0.44 


0.98 


0.34 


12.96 


0.13 


-1.69 


0.45 


-2.58 


0.54 


-2.06 


0.63 


16 


15.23 


1.99 


16.37 


1.18 


17.25 


1.02 


16.73 


0.84 


1.49 


0.98 


13.74 


0.34 


-1.13 


1.04 


-2.01 


1.17 


-1.50 


1.32 


17 


13.72 


1.26 


14.77 


0.77 


15.61 


0.67 


15.38 


0.56 


0.19 


0.08 


12.23 


0.17 


-1.05 


0.55 


-1.90 


0.64 


-1.66 


0.75 


18 


12.91 


1.63 


15.04 


1.02 


15.60 


0.88 


15.61 


0.73 


0.74 


0.42 


12.51 


0.16 


-2.13 


0.60 


-2.69 


0.73 


-2.70 


0.87 


19 


12.52 


2.45 


14.28 


1.52 


14.85 


1.32 


15.57 


1.10 


0.80 


0.68 


12.38 


0.20 


-1.76 


0.75 


-2.32 


0.91 


-3.04 


1.09 


20 


10.47 


0.55 


12.25 


0.34 


13.03 


0.29 


12.92 


0.24 


10.44 


1.99 


12.68 


0.06 


-1.77 


0.21 


-2.56 


0.25 


-2.45 


0.30 


21 


12.69 


1.24 


14.83 


0.77 


15.49 


0.67 


15.29 


0.56 


1.03 


0.45 


12.56 


0.14 


-2.14 


0.51 


-2.80 


0.61 


-2.60 


0.73 


22 


11.11 


6.43 


13.59 


4.00 


15.46 


3.46 


15.70 


2.93 


6.43 


14.42 


12.72 


0.77 


-2.48 


2.92 


-4.35 


3.59 


-4.59 


4.33 


23 


12.61 


1.05 


14.29 


0.65 


14.99 


0.57 


14.37 


0.47 


3.35 


1.22 


13.04 


0.13 


-1.68 


0.49 


-2.38 


0.60 


-1.76 


0.71 


24 


12.49 


0.86 


14.25 


0.54 


14.47 


0.46 


14.34 


0.39 


3.33 


1.00 


12.99 


0.10 


-1.77 


0.36 


-1.98 


0.43 


-1.86 


0.51 


25 


13.13 


0.68 


14.65 


0.41 


15.20 


0.35 


15.52 


0.31 


2.50 


0.58 


13.12 


0.15 


-1.52 


0.53 


-2.07 


0.63 


-2.39 


0.75 


26 


12.13 


1.56 


14.64 


0.97 


15.50 


0.84 


15.67 


0.70 


2.34 


1.28 


12.69 


0.14 


-2.51 


0.53 


-3.37 


0.64 


-3.54 


0.77 


27 


11.54 


1.83 


14.08 


1.14 


15.61 


0.99 


15.90 


0.82 


6.93 


4.43 


12.93 


0.17 


-2.54 


0.63 


-4.08 


0.78 


-4.37 


0.93 


28 


12.21 


0.51 


13.70 


0.32 


14.28 


0.28 


14.05 


0.23 


0.12 


0.05 


11.43 


0.19 


-1.49 


0.21 


-2.07 


0.25 


-1.84 


0.30 


29 


12.33 


1.08 


14.57 


0.67 


14.84 


0.58 


15.09 


0.49 


4.70 


1.77 


13.08 


0.10 


-2.24 


0.39 


-2.51 


0.47 


-2.76 


0.56 


30 


12.29 


1.40 


14.31 


0.87 


15.15 


0.76 


15.14 


0.63 


0.69 


0.34 


12.23 


0.13 


-2.03 


0.49 


-2.86 


0.60 


-2.85 


0.72 



■mOZ d^S 6 IA8ei6000/Md-ojjs^:ATXJ^: 



TABLE 3 — Continued 



Ap. 


152' 




U 




B 




V 




Lhc 


b,c 


log(Nt^, / 


L152)'' 


mi52 - 


mu 


mi52 - 


me 


mi52 - 1 


mv 




(mag.) 


i<^) 


(mag.) 


i<^) 


(mag.) 


i<^) 


(mag.) 


('^) 


(lum.) 


M 


(ratio) 


i'^) 


(mag) 


M 


(mag) 


i'^) 


(mag) 


M 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




(10) 




31 


12.51 


26.13 


16.21 


16.25 


17.31 


14.03 


19.18 


11.74 


2.59 


23.63 


12.89 


1.86 


-3.71 


7.09 


-4.80 


8.69 


-6.68 


10.54 


32 


11.82 


1.23 


14.01 


0.77 


14.78 


0.66 


15.07 


0.55 


1.88 


0.81 


12.48 


0.11 


-2.18 


0.41 


-2.96 


0.50 


-3.24 


0.60 


33 


13.36 


2.99 


15.94 


1.86 


16.91 


1.61 


17.15 


1.34 


1.26 


1.32 


12.92 


0.25 


-2.58 


0.94 


-3.54 


1.15 


-3.78 


1.37 


34 


10.39 


0.60 


12.85 


0.37 


13.77 


0.32 


13.94 


0.27 


11.19 


2.33 


12.68 


0.05 


-2.46 


0.19 


-3.38 


0.23 


-3.55 


0.28 


35 


12.18 


1.02 


14.39 


0.64 


15.09 


0.55 


15.19 


0.46 


2.15 


0.77 


12.67 


0.10 


-2.22 


0.37 


-2.91 


0.44 


-3.01 


0.52 


36 


11.24 


1.46 


13.65 


0.91 


14.73 


0.79 


14.80 


0.65 


9.70 


4.94 


12.96 


0.13 


-2.41 


0.48 


-3.49 


0.59 


-3.56 


0.71 


37 


12.22 


1.46 


14.42 


0.91 


15.58 


0.78 


15.78 


0.65 


6.16 


3.14 


13.15 


0.13 


-2.20 


0.49 


-3.36 


0.59 


-3.56 


0.71 


38 


11.44 


1.39 


14.04 


0.86 


14.89 


0.75 


15.46 


0.62 


3.61 


1.75 


12.60 


0.11 


-2.60 


0.41 


-3.45 


0.50 


-4.02 


0.60 


39 


12.57 


0.89 


14.75 


0.55 


15.44 


0.48 


15.37 


0.40 


3.64 


1.13 


13.06 


0.09 


-2.19 


0.33 


-2.87 


0.39 


-2.81 


0.46 


40 


13.24 


2.08 


20.72 


10.75 


16.43 


1.12 


16.76 


0.93 


0.57 


0.42 


12.53 


0.19 


-7.48 


10.69 


-3.19 


0.86 


-3.52 


1.03 


41 


10.59 


0.33 


12.86 


0.21 


14.02 


0.18 


14.04 


0.15 


20.35 


2.38 


13.02 


0.03 


-2.27 


0.12 


-3.43 


0.14 


-3.45 


0.17 


42 


12.93 


1.00 


14.64 


0.62 


15.76 


0.53 


15.79 


0.44 


5.10 


1.76 


13.35 


0.10 


-1.70 


0.36 


-2.82 


0.43 


-2.86 


0.51 


43 


14.22 


2.43 


15.42 


1.50 


16.51 


1.30 


16.71 


1.08 


1.53 


1.28 


13.34 


0.24 


-1.20 


0.81 


-2.29 


0.97 


-2.50 


1.13 


44 


11.66 


0.51 


13.70 


0.32 


14.64 


0.28 


14.76 


0.23 


7.32 


1.32 


13.00 


0.05 


-2.04 


0.17 


-2.98 


0.20 


-3.10 


0.24 


45 


13.74 


1.88 


16.08 


1.17 


16.79 


1.01 


16.90 


0.84 


1.60 


1.05 


13.17 


0.19 


-2.33 


0.70 


-3.05 


0.83 


-3.15 


0.98 


46 


12.83 


0.65 


14.71 


0.40 


15.89 


0.35 


15.80 


0.29 


3.83 


0.87 


13.19 


0.07 


-1.88 


0.25 


-3.06 


0.30 


-2.97 


0.35 


47 


10.55 


0.47 


12.75 


0.29 


13.94 


0.25 


14.06 


0.21 


14.52 


2.37 


12.86 


0.04 


-2.20 


0.16 


-3.38 


0.20 


-3.51 


0.24 


48 


12.88 


2.05 


15.51 


1.28 


16.20 


1.10 


16.87 


0.92 


1.28 


0.92 


12.73 


0.16 


-2.62 


0.60 


-3.32 


0.73 


-3.99 


0.87 



N4631'^ 



5.67 



0.01 



7.57 



0.01 



8.14 



0.01 



7.94 



0.00 



503.66 



12.44 



0.00 



-1.90 



0.01 



-2.47 



0.01 



-2.27 



0.01 



"mag. = -2.5 X log(f) -|- zoro.point, where zero points are U: -20.92, B: -20.45, V: -21.1 from Zombeck (1990) and 152: -21.1 from Stecher et al. (1992) 

Luminosities and Lyman continuum photon flux computed assuming distance to NGC 4631 is 7.5Mpc 
^Ha luminosity in units of 10 ergs see" 

Computed from flux in large aperture 



TABLE 5 
Derived Quantities 

















CSF 














Ap. 


IB Agc° 




Masses, 


IB Models'" 




M^^g^ 


Turn-ofF 




Star Formation Rates, 


CSF Models'' 




Ma^g"^ 




152 


U 


B 


V 


152 


U 


B 


V 


Ha 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


2.68 


0.08 


0.36 


0.52 


0.63 


0.40 


0.00 


964 


3160 


3550 


3730 


1070 


2490 


2 


4.23 


0.14 


0.19 


0.20 


0.24 


0.19 


2.24 


1880 


2760 


2680 


2940 


1860 


2430 


3 


4.52 


0.09 


0.06 


0.02 


0.01 


0.04 


2.46 


1040 


832 


264 


232 


1030 


679 


4 


5.43 


0.91 


1.92 


1.42 


1.33 


1.40 


3.25 


10700 


22100 


18800 


19100 


10600 


16300 


5 


5.81 


0.50 


0.92 


0.80 


0.80 


0.76 


3.65 


6360 


11300 


11000 


11600 


6470 


9360 


6 


10.19 


1.30 


1.64 


1.35 


1.17 


1.36 


8.20 


42000 


55200 


46600 


38600 


41400 


44800 


7 


3.06 


0.32 


0.66 


0.62 


0.66 


0.56 


1.00 


5480 


10400 


8400 


8030 


5510 


7560 


8 


4.14 


0.03 


0.05 


0.08 


0.07 


0.06 


2.19 


393 


659 


1130 


873 


390 


689 


9 


2.58 


0.06 


0.17 


0.28 


0.33 


0.21 


0.00 


675 


1360 


1790 


1760 


821 


1280 


10 


6.14 


0.48 


0.54 


0.48 


0.47 


0.49 


3.96 


6930 


7590 


5890 


5160 


6980 


6510 


11 


4.25 


0.07 


0.11 


0.08 


0.07 


0.08 


2.26 


878 


1550 


1050 


922 


869 


1050 


12 


6.82 


0.48 


0.52 


0.39 


0.39 


0.44 


4.55 


6910 


8710 


7150 


7070 


6880 


7340 


13 


4.18 


0.04 


0.03 


0.12 


0.06 


0.06 


2.21 


506 


486 


1620 


756 


502 


775 


14 


8.76 


2.41 


1.85 


1.44 


1.24 


1.74 


6.94 


58500 


46400 


36700 


29100 


58900 


45900 


15 


3.64 


0.05 


0.08 


0.06 


0.09 


0.07 


1.68 


593 


1410 


1290 


2170 


583 


1210 


16 


0.00 


0.02 


0.11 


0.14 


0.27 


0.13 


0.00 


108 


459 


436 


747 


415 


433 


17 


7.82 


0.14 


0.52 


0.43 


0.51 


0.40 


5.64 


2490 


8780 


7000 


8280 


2520 


5810 


18 


6.27 


0.21 


0.25 


0.27 


0.28 


0.25 


4.07 


2890 


4000 


4520 


4480 


2910 


3760 


19 


7.07 


0.36 


0.59 


0.58 


0.29 


0.45 


4.78 


5400 


10300 


11200 


5650 


5430 


7590 


20 


5.34 


1.65 


3.42 


2.83 


2.96 


2.72 


3.16 


19000 


38100 


36400 


41300 


18800 


30700 


21 


5.93 


0.24 


0.35 


0.32 


0.37 


0.32 


3.78 


3190 


4410 


4600 


5560 


3250 


4200 


22 


5.13 


0.86 


0.89 


0.27 


0.21 


0.56 


2.93 


9610 


10200 


3640 


3010 


9700 


7230 


23 


3.17 


0.09 


0.23 


0.29 


0.59 


0.30 


1.13 


1490 


3660 


4030 


7520 


1480 


3630 


24 


3.46 


0.13 


0.29 


0.46 


0.56 


0.36 


1.46 


1780 


3980 


6800 


8060 


1750 


4470 


25 


2.88 


0.05 


0.18 


0.28 


0.25 


0.19 


0.38 


801 


2340 


3010 


2380 


795 


1870 


26 


5.28 


0.35 


0.37 


0.29 


0.23 


0.31 


3.09 


3990 


4070 


3640 


3190 


4150 


3810 


27 


3.79 


0.33 


0.34 


0.18 


0.16 


0.25 


1.86 


4680 


4990 


2520 


2020 


4650 


3770 


28 


13.16 


1.22 


3.39 


3.94 


3.96 


3.13 


12.90 


_ 


_ 


_ 


_ 


_ 


_ 


29 


2.98 


0.11 


0.17 


0.33 


0.31 


0.23 


0.82 


1800 


2680 


4420 


3720 


1790 


2880 


30 


7.85 


0.53 


0.80 


0.68 


0.74 


0.69 


5.68 


9400 


13500 


10900 


10500 


9520 


10800 
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CSF 














Ap. 


IB Ago" 




Masses, 


IB Models'' 




M„„,,^ 


Turn-ofr= 




Star Format: 


ion Rates, 


CSF Models'* 




M„„/ 




152 


U 


B 


V 


152 


U 


B 


V 


Ha 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


31 


4.01 


0.15 


0.05 


0.05 


0.01 


0.06 


2.09 


2030 


730 


549 


102 


2020 


1090 


32 


6.48 


0.60 


0.68 


0.57 


0.44 


0.57 


4.25 


8510 


11100 


9960 


7480 


8430 


9090 


33 


3.84 


0.06 


0.06 


0.06 


0.05 


0.06 


1.92 


881 


907 


771 


648 


879 


817 


34 


5.35 


1.78 


1.96 


1.44 


1.16 


1.59 


3.17 


20500 


21900 


18500 


16200 


20500 


19500 


35 


5.37 


0.35 


0.48 


0.43 


0.37 


0.40 


3.19 


4020 


5360 


5540 


5180 


3930 


4810 


36 


3.63 


0.46 


0.42 


0.25 


0.23 


0.34 


1.67 


5860 


7170 


5510 


5410 


5760 


5940 


37 


2.80 


0.13 


0.26 


0.24 


0.24 


0.22 


0.08 


1760 


2780 


2040 


1810 


1760 


2030 


38 


5.72 


0.74 


0.69 


0.54 


0.31 


0.57 


3.56 


9300 


8450 


7350 


4370 


9330 


7760 


39 


3.06 


0.09 


0.15 


0.19 


0.24 


0.17 


0.99 


1500 


2330 


2610 


2930 


1510 


2180 


40 


6.15 


0.14 


0.00 


0.17 


0.14 


0.11 


3.97 


2050 


207 


2060 


1520 


2070 


1540 


41 


3.30 


0.67 


0.94 


0.62 


0.63 


0.71 


1.28 


9820 


13900 


10100 


10400 


9700 


10800 


42 


2.30 


0.09 


0.33 


0.32 


0.39 


0.28 


0.00 


902 


2260 


1720 


1770 


1410 


1610 


43 


2.33 


0.03 


0.16 


0.16 


0.16 


0.13 


0.00 


277 


1100 


859 


759 


424 


684 


44 


3.38 


0.25 


0.52 


0.56 


0.60 


0.48 


1.38 


3730 


6510 


5720 


5420 


3670 


5010 


45 


2.75 


0.03 


0.06 


0.09 


0.10 


0.07 


0.00 


427 


600 


663 


639 


444 


555 


46 


2.71 


0.08 


0.23 


0.21 


0.29 


0.20 


0.00 


989 


2120 


1520 


1760 


1060 


1490 


47 


4.24 


0.98 


1.29 


0.94 


0.94 


1.04 


2.25 


12800 


18400 


12700 


11800 


12700 


13700 


48 


5.08 


0.17 


0.14 


0.13 


0.07 


0.13 


2.88 


1850 


1720 


1820 


1010 


1870 


1650 



N4631'= 



20.03 



29.43 



26.10 



25.64 



25.30 



298000 



393000 



341000 



317000 



298000 



329000 



"Units: Myr 

''Units: lO'^ Mq 

"^Units: Myr, star formation initiated at 14 Myr 

''Units: Mg Myr"' 

'^Aperture Sum 



